Mechanical noise dependent Aging and Shear Banding behavior 
of a mesoscopic model of amorphous plasticity 
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We discuss aging and localization in a simple "Eshelby" mesoscopic model of amorphous plas- 
ticity. Plastic deformation is assumed to occur through a series of local reorganizations. Using a 
discretization of the mechanical fields on a discrete lattice, local reorganizations are modeled as 
local slip events. Local yield stresses are randomly distributed in space and invariant in time. Each 
plastic slip event induces a long-ranged elastic stress redistribution. Mimicking the effect of aging, 
we focus on the behavior of the model when the initial state is characterized by a distribution of 
high local yield stress values. A dramatic effect on the localization behavior is obtained: the system 
first spontaneously self-traps to form a shear band which then only slowly broadens. The higher the 
"age" parameter the more localized the plastic strain field. Two-time correlation computed on the 
stress field show a divergent correlation time with the age parameter. The amplitude of a local slip 
event (the prefactor of the Eshelby singularity) as compared to the yield stress distribution width 
acts here as a mechanical effective temperature- like parameter: the lower the slip increment, the 
higher the localization and the decorrelation time. 



I. INTRODUCTION 

While metals are characterized by a low elastic limit 
and large deformation before failure, their amorphous 
counterparts, Bulk Metallic Glasses (BMG) are known 
for the exact opposite, high mechanical strength and low 
ductility. The propensity of plastic deformation in BMG 
to localize and form shear bands is the main mechanism 
leading to mechanical failure. Understanding and thus 
controlling shear band formation is the main challenge 
that has so far limited the use of glasses as structural 
materials^ 

Although plastic deformation may be responsible for a 
significant heat production, recent experimental studies 
have shown that the origin of nucleation and propagation 
of shear bands could not be attributed to an adiabatic 
shear banding instability via local temperature rise^. In 
absence of such a thermal softening mechanism, Falk, 
Shi and collaborators have explored by Molecular Dy- 
namics simulations the hypothesis of a structural soften- 
ing mechanism 4 . They showed in particular that shear- 
banding was facilitated by a high degree of relaxation of 
the glassy structure. 

They proposed an interpretation of this phenomenon 
in the framework of the Shear Transformation Zones 
(STZ) theor^ describing plastic deformation of metal- 
lic glasses as resulting of local inelastic transformations 6 . 
Let us recall that according to STZ theory plas- 
tic deformation is assumed to result from a series 
of local reorganizations occurring within a population 
of "small" atomic/molecular clusters (zones) through 
micro-instabilities. Plastic deformation directly results 
from the balance between flips in the positive and neg- 
ative directions of these Shear Transformation Zones at 
rate that depends on an intensive parameter (e.g. free 
volume or effective structural temperature). Shi and Falk 



could associate associate the shear-band with a structural 
signature characterized by an effective temperature, re- 
flecting the higher potential energy in the band than 
in the still surrounding. In a similar spirit, Manning, 
Langer and collaborators^ proposed an enriched version 
of the STZ theory able to capture strain localization. The 
introduction of a relaxation-diffusion equation of the ef- 
fective structural temperature was in particular shown 
to induce shear-banding in aged structure (low effective 
structural temperature) and/or high shear-rate condi- 
tions. 

Independently, starting from the trap model devel- 
oped by Bouchaud 8 for the glass transition, Sollich, 
Cates and Lequeux^ developed a Soft Glassy Rheology 
(SGR) model to capture the rheology of complex flu- 
ids. In the trap model a landscape of traps of depth E 
drawn from an exponential distribution exp(— E/Eq) is 
assumed. A break-down of ergodicity naturally emerges 
at To = Eo/k. From this simplified view of the glass 
transition, Sollich et al introduce the mechanical stress 
as a bias to the energy landscape. It is important to note 
that the temperature in their model is not associated to 
a real thermal bath but is assumed to emerge from some 
mechanical noise a priori related to elastic interactions 
induced by local reorganizations. 

While STZ and SGR models capture part of the 
rich phenomenology of amorphous visco-plasticity, their 
mean-field character does not allow them to account 
for localization unless an additional ingredient is intro- 
duced. The latter can be the relaxation/diffusion of a 
state variable as discussed above and/or the inclusion of 
anisotropic elastic effect of local plastic events (Eshelby 
inclusion^) in the modeling. 

Building on the latter grounds several authors have 
developed "Eshelby" mesoscopic models to study plas- 
ticity of amorphous materials^ - — . Except in the case of 
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Ref J£ where a state variable is implemented or of Ref 
where the presence of walls traps plastic deformation, in 
such models, localization appears to be only transient 
and complex spatio-temporal correlations very similar to 
those observed in atomistic simulations emerge from the 
competition between diffusion and localizatio n 1 ^ 18 . 

Recently Fielding and collaborators^^ investigated 
an age-dependent transient shear banding behavior in 
different models where the shear banding was not trig- 
gered by an elastic or viscous softening constitutive law, 
but rather through an aging/rejuvenation behavior where 
the diffusive character of an internal variable would dic- 
tate the widening and progressive vanishing of an initial 
shear band. The introduction of such a mechanism in a 
variant of the SGR model results in a very slow ( "glassy" ) 
spreading of such shear bands. 

This age dependence of shear banding and its fast or 
glassy relaxation motivates us to reassess the question 
of the connection to be made between the glass theory 
inspired SGR model and the STZ model built from the 
identification of the microscopic mechanism of plastic- 
ity in amorphous materials. In particular, it has re- 
mained so far difficult to give a microscopic justification 
to the effective mechanical temperature defined in the 
SGR model 4 .* 2 !. 

In the following we present results about aging and 
localization obtained with the original mesoscopic model 
of plasticity presented in details in Ref We discuss in 
particular the effect of two parameters of the model which 
will appear to respectively mimic the age of the system 
before shearing and a mechanical effective temperature. 



II. DEFINITION OF THE MODEL 

Let us briefly recall the definition of the model (see 
Refill for more details). The mechanical fields are dis- 
cretized on a square lattice with a mesh size significantly 
larger than the typical scale of a plastic reorganization. 
Periodic boundary conditions are considered. The ma- 
terial is assumed to be elastically homogeneous, so that 
stresses and elastic moduli are scaled so that the steady- 
state local yield stress is unity. A local criterion of plas- 
ticity is considered. The initial distribution of local yield 
stress is denoted Pi(cr c ). Every time a local plastic crite- 
rion is satisfied at point ceo, a local slip Ae p occurs (we 
assume here that local plastic strains obey the same sym- 
metry as the external loading, pure shear in the present 
case, so that a simple scalar yield criterion can be cho- 
sen) with a random amplitude d drawn from a statistical 
distribution Q(d), Ae p (x) = dSo(x — Xo) where Sd is 
the Dirac distribution. Note that d is the product of 
the mean plastic strain by the "volume" of the trans- 
formation zone. This local slip d induces a long range 
redistribution of elastic stress with a quadrupolar sym- 
metry (see Ref j 17 i 22 for analytical and numerical details 
about this elastic propagator) Aa e i(x) = dG(x — Xo) 
with G(r, #) ~ Ad cos 40 jr 2 where A is the dimensionless 



elastic constant, r and the polar coordinates. The slip 
amplitude d, is drawn from a uniform distribution in the 
range [0; do]. 

After slip, the microstructure of the flipping zone has 
changed and a new value of the local yield stress is drawn 
from a distribution Ps(cr c ). The system is driven with 
an extremal dynamics so that only one site at a time is 
experiencing slip. The originality of the present depin- 
ning models relies in the anisotropic elastic interaction. 
Within this framework of dynamic phase transition, the 
choice of extremal dynamics ensures to drive the system 
at the verge of criticality: the macroscopic yield stress is 
given by the critical threshold. 

The above model may be seen as a depinning model for 
amorphous plasticity with a peculiar (anisotropic) elastic 
interaction. While the richness of the physics of depin- 
ning models mainly relies on the competition between 
elasticity and disorder, we see here that the anisotropic 
character and the abundance of soft modes in the elastic 
interaction which characterize the present model of amor- 
phous plasticity, naturally induce an additional competi- 
tion between localization and disorder. 

An implicit assumption used in our model is that the 
statistical distribution Ps(cr c ) used to renew the local 
plastic threshold under shear (i.e. after local slip) is 
the very same as the distribution of plastic thresholds in 
the initial configuration Pi(a c ). This hypothesis may be 
questioned. Indeed, various experimental and numerical 
results obtained in friction or in shearing granular mate- 
rial or complex fluids ^ 23 ! 24 seem to indicate an effect of 
the preparation of the material upon its behavior under 
shear. One may think for instance at the effect of den- 
sity of granular material: a loose (dense) packing tends 
to exhibit hardening (softening) while under shear their 
density progressively evolves toward a so-called "critical" 
value. 

In order to test the effect of our hypothesis we give in 
the following a bias to the initial thresholds distribution 
and try to test its consequences. Practically speaking, 
the yield stress distributions Pi(cr c ) (initial state) and 
Ps(&c) (under shear) are chosen as uniform in the ranges 
[5; 1 + S] and [0; 1] respectively. A positive (negative) 
value of S is expected to induce some softening (hard- 
ening) behavior since all threshold values above unity 
(below zero) should eventually be replaced by thresholds 
within the interval [0, 1]. 

We focus in the following discussion on the effect of 
these two parameters do and S. In the view developed 
by Sollich et al. 9 , the parameter do which gives the am- 
plitude of the mechanical noise induced by the elastic 
interactions may be thought as analogous to the effec- 
tive mechanical temperature x in the SGR model. How- 
ever, it is to be emphasized that this mechanical "noise" 
is strongly inhomogeneous in space and displays strong 
temporal and spatial correlations, absent from the SGR 
model. The second parameter S measures the shift be- 
tween the initial yield stress distribution, that uniform 
in [5; 1 + 5] and the distribution of new local yield stress 
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FIG. 1. Maps of plastic strain obtained from left to right at (e p ) — 1/16, 1/4, 1, 4 and 16 and from top to bottom with a bias 
value 5 — 0, 0.5 and 0.7 with a slip increment d — 0.3. 



under shear, uniform in the range [0; 1] may be related to 
the initial state of the system prior to shearing. Indeed 
it is expected that the older the glassy system, the more 
stable and the more difficult it is to shear. This effect 
is described here through a mere penalty in the initial 
yield stress. High mean values of the plastic thresholds 
should thus be associated with aged configurations of the 
glass. As discussed in 25 , a logarithmic increase of the 
yield stress with the age of the system is often observed 
in glassy materials: 5 « sq(T) log(t w /to) where to is a 
microscopic time scale. According to this perspective, 
the age of the system would simply be related to the bias 
6 through an exponential dependence. Yet another inter- 
pretation would consist of relating the parameter S to a 
structural temperature as discussed in Ref The two 

quantities are expected to vary in opposite ways when the 
structure relaxes. The more relaxed the glass, the higher 
5 and the lower the structural temperature. We now test 
these simple ideas against numerical simulations. 



III. MAPS OF PLASTIC ACTIVITY 

In Fig. [1] and [2] the evolution of the spatial dis- 
tribution of plastic strain under shear is shown for 
different values of the age-like parameter S (Fig. [1]) 
and of the temperature-like parameter do. We show 



snapshots of the plastic strain field taken at (e p ) = 
1/16, 1/4, 1, 4, 16. The value of the local strain is rep- 
resented with a grey scale, (the darker, the larger plastic 
strain). 

In Fig. [TJ the values do = 0.3 of the slip increment 
has been used. The first row corresponds to the value 
5 = 0. When using this un-aged initial configuration, 
we see that plastic strain first self-organizes along shear 
bands at ±7r/4 i.e., according to the maximum shear di- 
rections. This localization is however not persistent and 
after a transient, these shear bands diffuse throughout 
the system. The evolution obtained with a bias value 
8 = 0.5 (second row) is markedly different. Again plastic 
deformation first tends to form shear bands according to 
directions at ±7r/4, but remains essentially trapped in 
a strongly localized state. The formed shear band only 
slowly widens with "time" (mean plastic strain). The 
evolution obtained with a bias value S = 0.7 (third row) 
is very similar : formation of a persistent shear band be- 
fore an apparent diffusive widening of the band. Local- 
ization appears to be more intense and widening slower 
with this higher value of the age-like parameter S. 

Let us only note here that the way the plastic activity 
gets localized along a band is somewhat reminiscent of 
the behavior of an earlier model proposed by Torok and 
Roux 27 . In this study, the authors made evidence for a 
weak breaking of ergodicity which they relate to the pro- 



gressive building in the thresholds landscape of a valley 
(along the shear band) surrounded by ridges elevating 
significantly above the base level. Plastic activity thus 
tends to be confined in the valley and can no longer fully 
explore the disordered landscape. 

In Fig. [2j the values 5 = 0.5 of the age-like param- 
eter has been used. From top to left, the evolution of 
the plastic strain field is shown for values of the slip 
increment do = 0.03, 0.1, 0.3. A similar behavior as 
above is obtained. We see that the higher value of the 
temperature like parameter do, the less intense the local- 
ization and the faster the subsequent widening process. 
Age- and mechanical temperature-like parameters S and 
do thus seem to behave as could be expected, at least 



phenomenologically. 

IV. SLOW RELAXATION OF 
SHEAR-BANDING 

The residual stress field, is the self-balanced stress field 
which results from the local slip events taking place from 
the initial (stress free) state. The latter has a zero vol- 
ume average. It allows one (in conjunction with the local 
random yield threshold) to characterize the propensity 
of a site to undergo a plastic slip. This motivates the 
recourse to standard tools used for aging behavior char- 
acterization. Two-point correlation functions based on 
the residual stress field are proposed: 



^ ( \ _ {p'res {&w i %)&res %))x /-, \ 

{{(J res {s w , x)(7 res (£ w , x))x {p'res (%>) %)&res i^pi %)) x) 



where the symbol {...) x designates a spatial average over 
x. Note that the model does not depend on time as such; 
the global plastic strain plays the role of an evolution 
parameter. 

Such two-point correlation functions can be used to fol- 
low the formation and the subsequent relaxation of shear- 
banding^. In the following we only discuss the relaxation 
stage after the initial transient and full formation of the 
shear band. In Fig. [3] we present the dependence of the 
stress correlation for s w = 1; at this deformation level, 
which corresponds to the typical amplitude of the local 
plastic threshold, localization (if any) is fully set. 

The left panel of Fig. [3] shows the effect of the "age" 
parameter, with 5 = 0, 0.2, 0.4 ,0.6 ,0.8 and d = 0.2. 
In the un-aged configuration ((5 = 0), the system decorre- 
lates after a typical plastic strain s p = d^ 5 . This reflects 
the non persistence of localization in the standard un- 
aged case. In the case of an aged initial configuration 
we obtain significantly different results. The systems 
appears to decorrelate only after a plastic deformation 
growing exponentially with the parameter S. Moreover 
when fitting data with a simple stretched exponential, 
the exponent can be shown to transit from values slightly 
below unity in the un-aged case to values close or below 
1/2 in the more aged configurations. The shear banding 
persistence thus seems to directly depend on the age. 

Pursuing the above discussed analogy we now show in 
the center panel of Fig. [3] the correlation functions ob- 
tained with a fixed age parameter S = 0.6 for values of 
the slip increment parameter varying from do = 0.02 to 
do = 1 (computations were performed on lattices of size 
64 x 64 with 20 to 200 realizations). As could be antic- 
ipated from the above displayed maps of plastic defor- 
mation, the shear-banding persistence tends to increase 
inversely with the slip increment parameter, the lower 



do the higher the decorrelation time. The slip increment 
parameter do thus seems reasonably to act as the ampli- 
tude of a mechanical noise allowing the system to escape 
its trapped state. In other words, do which stands here 
for the product of the volume of a flipping zone times its 
typical plastic strain seems to be a good candidate for 
the elusive effective mechanical temperature discussed in 
the SGR modefi. 

We try to rationalize in the right panel of Fig. [3]the age 
and mechanical noise dependence of the shear-banding 
persistence. Exploring the two-dimensional space of pa- 
rameters S and do, using a simple stretched exponen- 
tial fitting procedure, we extracted the typical plastic 
strain £* associated with stress decorrelation after shear- 
banding formation (e w = 1 in the above notations). This 
allows us to propose a reasonable scaling dependence: 

where 

cp(x -> 0) « A and cp(x -> oo) « C.e Bx . 

The choice a = 0.5 and b = 0.2 allowed us 
to obtain a reasonable collapse of the data col- 
lected for do e [0.02,0.03,0.05,0.1,0.2,0.5,1] and 5 £ 
[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9]. In Fig. [3] an in- 
dicative exponential curve is shown to account for the 
high age and/or low mechanical noise shear banding slow 
relaxation behavior. 



V. CONCLUSION 

To summarize, we showed that our simple Eshelby- 
like mesoscopic model of amorphous plasticity exhibits a 
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FIG. 2. Maps of plastic strain obtained from left to right at (s p ) — 1/16, 1/4, 1, 4 and 16 with a bias value 5 — 0.5, and from 
top to bottom with a slip increment do = 0.3, 0.1 and 0.03. 
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FIG. 3. Effect of the age-like parameter 5 (Left) and of the mechanical noise amplitude do (Center) on the two-point stress 
correlation function with e w = 1. (Right) Dependence on S and do of the typical plastic deformation e* needed to relax shear 
banding. After rescaling data can be reasonably collapsed onto a single master curve; the dashed line indicates an exponential 
behavior accounting for the glassiness of shear-banding at high 8 - low do. 



striking dependence on initial conditions. The introduc- 
tion of a simple bias to shift the initial distribution of lo- 
cal yield stress values from its counterpart used to renew 
the yield stress after local reorganization has a remark- 
able consequence: the systems self-traps in a localized 
state to form a shear band and remains so for a longer 
and longer "time" when the bias value increases. This 
bias can thus be interpreted as an estimator of the age of 
the system before shearing (or related to some effective 
structural temperature^. 

Moreover we show that the ratio of the typical slip 



increment (more rigorously in the formalism of the Es- 
helby inclusion, the volume of a reorganizing zone times 
its typical plastic strain) on the typical plastic yield stress 
acts as an effective mechanical temperature in the sense 
proposed in the SGR model of Sollich et alA This pa- 
rameter indeed gives the amplitude of the mechanical 
noise induced by successive reorganizations. The lower 
this amplitude, the longer the systems gets trapped and 
the slower the widening of the shear bands. 

In conclusion the present depinning model of amor- 
phous plasticity appears to reproduce shear-banding, a 
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crucial feature of the phenomenology of metallic glasses. 
The nucleation step of the shear-band is followed by a 
slow broadening step of the band. The latter is quanti- 
tatively characterized by a slow relaxation of the stress- 
stress correlation. Note that this slow dynamics spon- 
taneously emerges in absence of any prescribed internal 
relaxation timescale. 

The behavior of the model is controlled by two pa- 
rameters that can be associated to a structural effective 
temperature^) and a mechanical effective temperature^ 



respectively. This model may thus contribute to clarify 
the respective effects of structural relaxation and me- 
chanical noise induced by local reorganizations in the 
plastic behavior of amorphous materials. 
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